Magnetic excitations in Kondo liquid: Superconductivity and Hidden Magnetic 

Quantum Critical Fluctuations 
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We report Knight shift experiments on the superconducting heavy electron material CeCoIns 
that allow one to track with some precision the behavior of the heavy electron Kondo liquid in the 
superconducting state with results in agreement with BCS theory. An analysis of the ^^^In nuclear 
quadrupole resonance (NQR) spin-lattice relaxation rate T^^ measurements under pressure reveals 
the presence of 2d magnetic quantum critical fluctuations in the heavy electron component that are 
a promising candidate for the pairing mechanism in this material. Our results are consistent with 
an antiferromagnetic quantum critical point (QCP) located at slightly negative pressure in CeCoIns 
and provide additional evidence for significant similarities between the heavy electron materials and 
the high Tc cuprates. 

PACS numbers: 74.70.Tx, 76.60.-k, 76.60.Cq 



CD 



O 

o 



> ■ 



o 

On 
O 



- '—I 

X 
J3 



A central issue in heavy electron physics is the origin 
of unconventional superconductivity (SC) at low temper- 
atures, and the role played by quantum critical fluctua- 
tions in bringing that about. Valuable information on 
both is provided by nuclear magnetic resonance (NMR) 
experiments of the Knight shift and spin-lattice relax- 
ation rates but as is so often the case the devil is 
in the details - in this case being able to separate out 
the contributions made by local moments to both quan- 
tities, so that one can be confident one is measuring the 
behavior of the heavy electrons alone and not a combi- 
nation of the two. We are now in a position to do this, 
thanks to a phenomenological two-fluid description that 
treats the emergence of heavy electrons as the formation 
of a distinct new quantum state of matter that displays 
universal behavior below a characteristic temperature T* 
0, & i, i, H- The growth of this new phase is accom- 
panied by a loss in strength of the spectral weight of 
the local moment contribution, and both phenomena are 
characterized by an order parameter, /(T). so that it be- 
comes possible, for a given physical phenomenon, to track 
the role played by residual local moments and focus on 
the behavior of the emergent phase, the heavy electron 
Kondo liquid. 

In this communication, we first apply the two-fluid 
description to new measurements of the Knight shift in 
CeCoIns, and show that it provides a natural explanation 
for the anomalous behavior of the NMR Knight shift that 
is seen in both the normal and superconducting phases. 
The Knight shift anomaly, defined as the deviation of 
the total Knight shift from the magnetic susceptibility, 
shows scaling behavior in both phases, is independent of 
the probe nuclei, follows the predicted Kondo liquid uni- 
versal behavior in the normal phase [li, 13, and BCS 
theory in the superconducting phase [7[, thereby estab- 
lishing unambiguously that the heavy electron SC origi- 



nates in the condensation of the heavy Kondo liquid. 

We then carry out a two-fiuid analysis of the NQR mea- 
surements of the ^"'^^In spin-lattice relaxation rate in the 
normal state of CeCoIns Q, and find that for pressures 
up to 1.2 GPa, the heavy electron spin-lattice relaxation 
rate, T^j^ , takes the simple form expected for 2d magnetic 
quantum critical fiuctuations, Ti^T oz T + Tq This 
behavior persists up to T*, while the offset, Tq, tends 
toward zero as one moves toward a candidate QCP at a 
slightly negative pressure [lo| . Such behavior is hidden 
in the overall spin-lattice relaxation rate by the presence 
of a substantial contribution from the unhybridized local 
moments, and can only be identified once this component 
is subtracted. 

We have performed NMR measurements on CeCoIns 
down to 0.04 K and the results in Fig. [1] complement our 
previous measurements According to the two-fluid 
scenario 0, 01 , the magnetic susceptibility and Knight 
shift are related through 



X = Xh+Xh 
K-Ko = Axh + Bxi, 



(1) 



where Xh is the heavy electron susceptibility of the Kondo 
liquid and xi is the contribution to the susceptibility from 
the residual unhybridized local moments. A and B are 
the hyperfine couplings of the probe nucleus to the two 
components, respectively. Kq is a constant offset that 
accounts for other contributions such as the Pauli term 
from the "light" conduction electrons. Above T* ^ 40 K, 
Xh — and the Knight shifts measured for the In(l); Co, 
and In(2) sites arc proportional to the magnetic suscep- 
tibility, K = Kq + Bx- Below T* . the emergence of Xh 
leads to the deviation from this linear relation shown in 
Fig- [Ha). The Knight shift anomaly is defined as 



= K-Ko-Bx = iA-B)xh, 



(2) 



and provides a direct measurement of Xh ■ An analysis of 
the Knight shift anomaly in many heavy electron mate- 
rials has led to the conclusion that Xh displays universal 
behavior In Fig.H the planar Knight shift anoma- 

lies at these probe nuclei are found to be in good agree- 
ment with the predicted universal Kondo liquid density 
of states i, phiT) = (1 - T/T*)^^^[l + In(rVr)]. 

We consider next the extension of the two-fluid sce- 
nario to the superconducting phase. The planar Knight 
shift of CeCoIns has several unexpected features that are 
not present in most other heavy electron materials and 
find a simple explanation within the two-fluid scenario. 
First, no anomaly is observed in the planar Knight shift 
data at the In(l) site. According to Eq. ([2]), this sug- 
gests a cancellation of the In(l) hyperfine couplings to 
the residual local moments and the heavy Kondo liquid. 
Since the superconducting transition does not change the 
hyperfine couplings, this cancellation is expected to hold 
across Tc- Therefore, the In(l) nucleus probes the to- 
tal spin susceptibility in the whole temperature range. 
This allows us to determine the intrinsic spin suscepti- 
bility, X = {K - i^o)in(i)/5in(i), even below Tc where 
the measured susceptibility goes negative. So a Knight 
shift anomaly can still be defined in the superconducting 
phase for probe nuclei other than In(l) by using Eq. ([2]). 
In general cases where A ^ B for all probe nuclei, the 
Knight shift anomaly or Xh can be obtained (up to a con- 
stant prcfactor) in a similar way even though the effective 
susceptibility defined as {K — Ko)/B docs not give the 
true spin susceptibility x- 

Second, the planar Knight shift at the In(2)^ site is a 
constant above T*, but evolves with temperature below 
T*. This indicates that the In(2)j^ nucleus is not coupled 
to the localized moments. This special feature allows us 
to identify unambiguously the Kondo liquid behavior in 
both the normal and superconducting phases and pro- 
vides an independent check on the Knight shift anomaly 
measured by other probe nuclei (In(2)|[ or Co). As can 
be seen in Fig.[2l the In(2)^ Knight shift anomaly above 
Tc does fall upon the universal curve for the Kondo liquid 
after performing a simple scaling. Without any artificial 
manipulation, the Knight shift data at In(2)j^ nucleus 
provides added confirmation for the existence of the pre- 
dicted universal Kondo liquid state. 

Fig. [IJb) compares our new experimental data in the 
superconducting phase with that obtained earlier 



By using the effective susceptibility obtained from the 
Knight shift at In(l) nucleus, the scaled Knight shift 
at In(2)j^, and the subtracted Knight shift anomalies at 
In(2)|| and Co are compared in Fig. [2] for the whole tem- 
perature range below T*. Due to the difference in the ap- 
plied external magnetic field, the superconducting tran- 
sition temperatures are slightly different. However, the 
similar scaling behavior at different probe nuclei supports 
the Kondo liquid scenario in both the normal and super- 
conducting phases and shows distinctly the condensation 
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FIG. 1: (Color online) Planar Knight shift of CeCoIns at dif- 
ferent probe nuclei in both normal (a) and superconducting 
(b) phases. The solid lines are the measured magnetic sus- 
ceptibility. The solid symbols denote previous experimental 
data from Ref. [lH and the open symbols are our new data. 
The difference in Tc is due to the different external magnetic 
field for the measurements. The inset indicates the position 
of the different probe nuclei (In(l), Co, In(2)x and In(2)||) in 
the unit cell. There are two inequivalent In(2) sites, depend- 
ing on whether the field is parallel (In(2)||) or perpendicular 
(In(2)x) to the unit cell face. 



of the Kondo liquid into the superconducting state. In- 
deed, the rapid decrease of the anomalies follow exactly 
the BCS prediction for d-wave superconductivity 0], 



(T) ~ KanorniO) K / dE 



df{E) 
dE 



N{E), (3) 



where f{E) is the Fermi distribution function and 
N{E) oc {\E\/^E^ - ^k{Tf)FS is the BCS density of 
states. The superconducting gap function is described by 
the usual interpolation formula 



Afc(r) =gfcA(0)tanh 



1 



Te (Tc 

dT ^ A(0)2 T 

^(4) 

where A(0) is the gap amplitude and describes the 
gap symmetry. In Fig. [2l the best fit for d-wave gives 
A(0)/rc ~ 4.5 and a specific heat jump AC/C ~ 4, both 
in good agreement with previous estimates [1, ■ We 
therefore conclude that the Kondo liquid is responsible 
for the heavy electron superconductivity. 
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FIG. 2: (Color online) Planar Knight shift anomaly of 
CeCoIns at different probe nuclei in both normal and super- 
conducting phases. The solid lines are fit to the universal 
Kondo liquid density of states in the normal phase Q and 
BCS theory in the superconducting phase 0]. The best fit 
gives a superconducting gap A{0) /Tc ~ 4.5 and a specific 
heat jump AC/C ~ 4, both in good agreement with previous 
estimates [lo|. 



The extension to the superconducting state reveals 
a third unexpected feature, a negative Knight shift 
anomaly at zero temperature, which is best understood in 
the In(2)j_ data. Above T*, In(2)j^ probes only conduc- 
tion electrons that have yet to hybridize. These electrons 
contribute a constant term in Kq in the normal phase. 
Therefore, we surmise that the negative anomaly relative 
to its high temperature offset Kq originates in the loss of 
these conduction electrons, in addition to the heavy elec- 
tron Xh, to the superconducting condensate. A similar 
loss of their contribution is expected for the other probes, 
In(2)|| and Co, and its presence in Fig. [2] confirms this 
physical picture. Since the In(2)|| and Co Knight shift 
anomalies are obtained after subtracting an effective sus- 
ceptibility deduced from the In(l) data, this picture also 
explains the upturn of the In(2)|| and Co Knight shifts 
below Tc and their consequential intersection with the 
In(l) Knight shift seen in Fig.[T]^b). 

Important information on the spin fluctuation excita- 
tions that may be responsible for the measured d-wave 
pairing may be obtained from the measured ^^^In NQR 
spin-lattice relaxation rate T^^ of CeCoIns under pres- 
sure Q shown in Fig. [31 Since, however, below T* both 
local moment and heavy electron spin fluctuations con- 
tribute to the overall relaxation rate, it is essential to 
separate these out before drawing any conclusion about 
the nature of either spin fluctuation spectrum. The two- 
fluid description makes this possible. Above T*, T-f^ 
measures the strength of the magnetic fluctuations of the 




FIG. 3: "^In NQR relaxation rate of CeCoIns at different 
pressures analyzed by the two-fluid scenario. The data are 
reproduced from [^. The solid lines are the best fit to the 
total Tj^^ and the dotted lines show the local moment contri- 
bution. The inset compares the magnitude of Ti at T* with 
the exchange interaction T* . 

unhybridized local moments while below T*, the local 
moment contribution is reduced by their hybridization 
and a heavy electron component emerges, producing the 
change in the temperature dependence of the total T^^ 
that has been seen in many other heavy electron materi- 
als Assuming no interference between the local and 
heavy electron contributions to Tj~^, the two-fluid model 
suggests that it takes the following form: 

1_ ^ 1 - f{T) fiTl 

where /(T) = (1 — T jT*)^!"^ is the order parameter and 
represents the fraction of the local moments that are con- 
verted into the itinerant Kondo liquid. The local term 
T\i is from the unhybridized local moments and its tem- 
perature dependence can be obtained by fitting to the 
high temperature T\ above T* and making the assump- 
tion that this temperature dependence persists below T* . 
The basis for our analysis is shown in Fig. [3] where the 
experimental data take the form T^^ = a + bT above T*, 
that has been seen for many heavy electron materials. It 
originates in the weak coupling between local moments 
that gives rise to their measured Curie- Weiss susceptibil- 
ity, and its magnitude near T* scales with the strength 
of the nearest neighbor inter-ion exchange interaction, 
which is of order T* Q . 

For CeCoIns at different pressures, we estimate T* 
from transport experiments displaying a Hall anomaly 
we find T* ~ 50 K at ambient pressure, ~ 60K at 
0.58 GPa and ~ 70 K at 1.2 GPa. We note that the NQR 
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FIG. 4: Kondo liquid TihT for CeCoIns at diflterent pres- 
sures. The linear temperature dependence of TihT signals 
the presence of quantum critical fluctuations. The inset plots 
the deduced To that measures the distance from an antiferro- 
magnetic quantum critical point. 



T[' contains the efTcct of both c and planar directions, so 
we should use an average T* that is slightly larger than 
the planar one used for the planar Knight shift analysis 
due to anisotropy. Assuming no change in the hyper- 
fine coupling with pressure, the inset of Fig. [3] compares 
Ti(T*) and T* and confirms their proportionality. 

To obtain the heavy electron contribution to the re- 
laxation rate, we assume that the intrinsic temperature 
dependence of the local moment contribution to T^i^ ex- 
tends to lower temperatures; what changes is its strength. 
The dotted lines in Fig. [3] give the local moment contri- 
bution [1 — f{T)]/Tii. The Kondo liquid contribution can 
then be obtained by using Eq. ([5]). The resulting TihT 
are plotted in Fig. [4l Interestingly, this Kondo liquid re- 
laxation rate fits well to the following simple formula for 
all pressures below 1.2 GPa, 



TihT cx (T + To) 



(6) 



where Tq is a positive constant of the order of a few 
Kelvin. The slope increases with increasing pressure or 
T* . If we take the above linear formula for the Kondo 
liquid as our starting point and fit to the experimental 
data, we do find a good fit (sohd lines in Fig. [3]). The scal- 
ing behavior in Eq. ^ has been found in many high Tc 
cuprates where it signals the presence of quantum critical 
fiuctuations of a nearly 2d spin liquid This provides 
additional evidence of the similarity between these two 
kinds of unconventional superconductors. The deviation 
from Eq. ^ above 40 K arises either from experimental 
errors or the breakdown of the quantum critical behavior 
in that high temperature region. 



To measures the distance from a quantum critical 
point. In the inset of Fig. |4l we plot Tq as a function of 
pressure and conclude that the quantum critical point for 
this material is at a slightly neg ative pressure, consistent 
with previous expectations [10| . The increasing strength 
of these quantum critical fluctuations f{T)/Tih of the 
Kondo liquid with lowering temperature for Tq < T < T* 
is in strong contrast to the overall decrease of Tj~^, while 
these are a highly promising candidate for their super- 
conductivity below Tc- 

To summarize, we find scaling behavior for the Knight 
shift anomaly in both the normal and superconducting 
phases of the heavy electron superconductor CeCoIns 
that provides a natural explanation for the various un- 
usual features of the Knight shift data and refiects di- 
rectly the heavy electron Kondo liquid condensation into 
the superconducting state. The quantum critical fluctu- 
ations we deduce from our two-fluid analysis of normal 
state behavior provide a highly promising candidate for 
the physical origin of their superconductivity, while the 
pressure variation of the calculated temperature depen- 
dence of the intrinsic heavy electron Kondo liquid spin- 
lattice relaxation rate suggests that the magnetic QCP 
is located in this material at a slightly negative pressure, 
as suggested by previous experiments. The theoretical 
method applied in this analysis can be easily extended 
to other heavy electron materials, and we find that on 
doing so for CeRhIn,5 that the magnetic quantum critical 
fluctuations seen there are likewise in the heavy electron 
channel, suggesting that its long range magnetic order is 
triggered by a spin density wave instability in this com- 
ponent. 
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